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, Short length quantum wires (quantum contacts) exhibit a conductance structure at the 

' value of conductance close to 0.7x2e^//i. The structure is also called the conductance anomaly. 

\ In longer contacts the structure evolves to the lower values of conductance. We demonstrate 

that this structure is related to the development of charge density waves within the con- 
^ tact. This is a precursor for Wigner crystallization. Many-body Hartree-Fock calculations of 

, conductance are performed. The results are in agreement with experimental data. 

• PACS: 73.61.-r, 73.23.Ad, 71.45.Lr 

(N 

, The quantized conducjbance G ~ nG2, n — 1,2, 3, G2 — 2e^ /h, through a narrow quantum point contact 

' was discovered in 1988lIi3. This quantization can be understood within a one- dimensional (ID) non-interacting 
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I. INTRODUCTION 



electron gas picture, see e.g. Ref.lj. In the present work we are interested in a deviation from thjCjnteger 
quantization. This deviation, the so called "0.7 anomaly" has been found in experimental worksu'u. The 
anomaly is a shoulder- like feature or a narrow plateau at G w O.7G2. The recent workQ demonstrates that 
the anomaly evolves down to G ~ O.5G2 in longer quantum contacts- Another highly interesting recent 
experimental development is an indication of the above barrier excitationQ. The excitation is probably related 

■ to the "0.7 anomaly". 

^ ■ There is practically no doubt that the structure is due to the interaction between the electrons. There 
O have been numerous attempts to explain the pQ-T anomaly" by spontaneous magnetization of the ID quantum 
1^1 ' wire due to exchange Coulomb interactionrtll. A somewhat similar, but still differpitjdea is to explain the 
anomaly by formation of a two-electron bound state with total spin 5=1, see Refs.Ej'tj. 

It is not clear why_the spontaneous magnetization scenario can be valid since it contradicts to the rigorous 
^ Lieb-Mattis theorenJii. The theorem claims that any ID many-body system with quadratic kinetic energy 
and with a potential interaction has the ground state with zero spin. The same argument is applicable to the 
S=l bound state picture. Strictly speaking in the later case one can avoid a formal contradiction with the 
, theorem saying that the bound state is localized somewhere in the transition region between the 2D or 3D 

■ conductor and the ID wire. However this scenario seems unlikely and in the present work I consider only S=0 
T-H [ ground state in accordance with Lieb-Mattis theorem. 

. The effect of Wigner crystallization in low density 2D and 3D electron gas is well known. There is no such 
> ■ a thing as an ideal ID Wigner crystal because quantum fluctuations destroy any long range order. However 
it was pointed out some time agotJ that impurities suppress the fluctuations and pin the ID Wigner crystal. 
Surprisingly this happens at pretty high density of electrons n: n k, 0.5/aB, where ub is the Bohr radius. 
A quantum contact (=short quantum wire) is not uniform and hence it pins the crystal. This must result 
(— ( in nontrivial behavior of the contact conductance. It is the effect we study in the present work using the 
Q Hartree-Fock (HF) method. We demonstrate that the effect results in the conductance anomaly at G ~ O.7G2 
O for a short contact and in several anomalies for longer contacts. In the present work we consider only zero 
temperature case. However in conclusions we comment on temperature dependence. 

The idea of relation between the conductance anomaly and the Wigner crystallization was also suggested 
in a recent paperil3. Authors of this paper discuss a ferromagnetic state embedded in the Wigner crystal. 
We follow the different way considering the ground state with total spin zero in accordance with Lieb-Mattis 
theorem. 

In independent particle approximation, i.e. in the case of an ideal electron gas, a calculation of the conduc- 
tance for a given transverse channel is straightforward. Assuming adiabatic transmission one gets 

2e2 

where T is the barrier transmission probability at Fermi energ In the case of interacting particles this 

formula is also valid because before and after the potential barrier the density of electrons is high enough. 
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and hence the interaction is neghgible. However one can not use the single particle description to calculate 
the transmission probability T because in the vicinity of the barrier the electron density is low, and hence the 
many-body effects are very important. The main question is how to calculate the transmission probability Tl 
To do this we apply the following trick. Consider the liquid of electrons on a ID ring with a potential barrier 
somewhere on the ring. There is no a current in the ground state of the system. Now let us apply a magnetic 
flux through the ring. This flux induces the electric current. On the one hand the current is related to the 
barrier transmission probability T. On the other hand, the current can be calculated using the HF method. 
This allows us to find T with account of many-body effects. 

Structure of the present paper is the following: The relation between the current on the ring and the trans- 
mission probability is derived in Sec. II. In Sec. Ill we demonstrate how the Hartree-Fock method describes 
development of the charge-density wave at low electron density. The results of selfconsistent calculations of 
the transmission probabilities for different barriers are presented in Sec. IV. Sec. V summarizes the results 
and presents our conclusions. 



II. RELATION BETWEEN THE BARRIER TRANSMISSION PROBABILITY AND THE 
CURRENT IN THE RING INDUCED BY THE MAGNETIC FLUX. 

Let us consider a single particle in a potential U{x). The potential is localized near x = 0, U{x) = at 
|a;| > 1/2. For simplicity we assume that the potential is symmetric, U{x) = U{—x). In the scattering problem 
the wave function outside the region of the potential is of the form 

x<-l/2: X = e''^''' + ^e-*'=^ (2) 
x>l/2: x = Be'''\ 

where k = \J2E is the momentum of the particle, A is the reflection amplitude and B is the transmission 
amplitude. Throughout the paper we set the electron mass, electric charge, and the Planck constant equal to 
unity, m — e = 'h=\. The transmission probability is T ~ One can also use the basis of standing waves. 

Since the potential is symmetric, there is one symmetric and one antisymmetric solution at a given energy E 

x < —1/2 : x± — cos{kx + S±), (3) 
X > 1/2 : x± = ± cos{kx — S±). 

Comparing eqs. (|^) and (Q) one can express the transmission and reflection amplitudes in terms of the 
scattering phases 

A = e"'^ cos (5, 

B = -ie-'^shiS, (4) 
T = \B\^ = sin^J, 



where 



Now we switch the vector potential A on 



A = 6+ + 6., (5) 



H=^^^ + U{x), (6) 

and close the ring imposing the periodic boundary condition 

ibi~L/2)^iP{L/2). (7) 

We will assume that the length of the ring L is much larger than the size of the region where the potential 
U{x) is nonzero: L ^ I. Note that the vector potential we introduce is the pure gauge one. We need it to 
induce the orbital current in the ring. The magnetic field is zero. The induced current depends only on the 
Bohm-Aharonov phase </> = AL. In the many-body problem the current has period = tt. The current is 
maximum at = ±7r/2, ±37r/2, .... For further considerations we take the first maximum 
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A = t:/2L, = 7r/2. 



(8) 



Solution of the Schrodinger equation with Hamiltonian ^ , and the periodic boundary condition (0) gives 
the following single electron wave functions outside of the region of the potential. 



v/i(r 



Here l/2<x < L- 1/2, 



a± 



1 =F *e 



1 ± ie-"-^ ' 

and the momentum k obeys the following quantization condition 



ikL I ■ iA 

e — ±te 



(9) 



(10) 



(11) 



There are two solutions: ^+ and To elucidate meaning of these solutions let us look at two limiting cases. 
The first limit is an infinitely high potential barrier, U oo. In this case the scattering phases are equal, 
5+ = S^. Substituting these phases into eqs. ([lO| ), one finds that the wave functions ^± coincide with 
symmetric and antisymmetric standing waves In the second limit there is no barrier, U — 0. Then the 
scattering phases are the following: S+ — 0, S- — —Tr/2. One can prove that in this case the states (||) describe 
clockwise and anti-clockwise rotations, ip-^- — > exp(ikx), — > exp(— iA;a; + 2iAx). 
The quantization condition (^l]) gives the following momenta 



= ^(m ± 1/4 + A/2^), m = 0, 1, 2, 3..., 



(12) 



and energy levels 



1 

2 
2 



2ti 



2ti 



(m + A/27r)^ , 



(13) 



— ) (to - 1/2 + A/27r)' 



Using the wave functions (H) one easily finds the electric current in each single particle state 

27r 



^ ra 



L2 



(to + A/27r) sin 5, 



(14) 



27r 



2 (to, - 1/2 + A/27r)sin5. 



Up to now we were considering the single particle problem. Now let us consider the many body problem. To 
do this we fill all the single particle states up to the Fermi energy. Each orbital state is filled by two electrons, 
spin up and spin down. Electric currents Jm'' and Jm ^ partly compensate each other, so the total electric 
current of the many-body system is 



J — 2 ^ ^ f + J^j 



An 

l2 



27r 

L2 



E 



dsin^ 
' dk 



^ (m + A/27r)((sin<5(fc(+)) - sin<5(fc(-))) + 1 sin<5(fc( 

d{k sin S) 



2tt sr-^ d{k sin S) 

l2 ^ dk 



dk 



dk — —— sm Of — -zr^ sm d f ■ 
L 2L^ 



(15) 



Here N is the number of electrons on the ring, kF — 7riV/2L is the Fermi momentum, and is the scattering 
phase (5, see eq. (|5|), taken at the Fermi energy. 

Eq. ( [TBI ) gives the main tool we use in the present work. First we solve numerically the many body problem 
on the ring without the barrier, [/ = 0, and find the electric current Jq- Second we solve numerically the many 
body problem with the barrier, U 0, and find the electric current Jij. Then using eqs. ( p^ ) and (^) we 
conclude that the transmission probability at the Fermi energy is 
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T = {Ju/Jof 



(16) 



This is exactly the transmission probabihty wc need for the conductance (|l|). 

In the derivation of eqs. ( [Tsl ) and ( [T^ ) we neglect electron-electron interaction outside the potential barrier. 
This is valid if the density of electrons is high enough. In the vicinity of the barrier the electron density is 
low and hence the electron-electron interaction is very important. This does not contradict to the presented 
derivation. The information about the electron-electron interaction is hidden in the scattering phases. 

Finally we would like to comment on the Luttinger liquid behavior that is related to the long-range fluc- 
tuations in the ID interacting Fermi gas. The HF method we use to analyze the electron ring certainly does 
not take into account the long-range fluctuations. However this is irrelevant to the calculation of the barrier 
transmission probability. The matter is that the number of electrons above the quantum barrier is not more 
than several, so there is no room for the long-range fluctuations within the barrier. The long ring outside the 
barrier is just a technical trick to describe the reservoir of electrons. There is no need to take into account the 
long-range fluctuations in the reservoir. 



III. HARTREE-FOCK (HF) METHOD FOR INTERACTING ELECTRONS ON THE RING IN 
THE EXTERNAL GAUGE FIELD. CHARGE-DENSITY WAVE PINNED BY THE IMPURITY. 



The Hamiltonian of the many body system we consider is of the form 

E [^^^-^ + Uix,)] + 1 E ^(^- ^.)' (17) 



H 
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where indexes i and j numerates electrons, Xi is the periodic coordinate {0 < x < L), A is the gauge field (^, 
U{x) is the potential barrier, and Vij is the electron-electron interaction. We use atomic units, so coordinates 
are measured in unites of Bohr radius, as = th^ /me^, and energies are measured in units of Eunit — "^^^/^rf^n 
where m is the effective electron mass and e is the dielectric constant. For experimental conditions of worksoQ 
these values are the following: as ~ 10~^/im, Eunit ~ \Q~^eV . The electron-electron Coulomb repulsion we 
take in the form 

V{x,y)^^= \ (18) 
^al + D^{x,y) 

where at is the effective width of the transverse-channel, and D{x^ y) is the length of the shortest arc between 
the points x and y on the ring, see commentta. The transverse motion of electrons is frozen due to some 
confining potential that can be approximated as 

U.ip) = (19) 

According to the dataBi the energy splitting between the transverse channels is Tiwj^ ~ 4 — 7meV. This gives 
the following value of the first channel transverse size 

at = l/^/rnu^l Ki 2aB- (20) 

Averaging over the oscillator wave function of the transverse motion shows that the effective ID Coulomb 
interaction is well approximated by eq. (p^. We stress that the_Yalue of at is close to the critical distance 
between the electrons at which the Wigner crystallization occursll3. This is a fortunate coincidence. Larger 
value of at would suppress the conductance anomaly, and this really happens in higher transverse channels. 
Substantially smaller value of at would require a theoretical technique more sophisticated than the HF method 
employed in the present work. We will discuss this point later. Let us also comment on the long-range behavior 
of the interaction (p^). Because of the ring geometry the interaction (|8|) at |x — y| ~ L is somewhat different 
from the Coulomb one. However the so long range tail of the interaction does not influence the results. Anyway 
in a real system there is a screening by gates and many other effects that influence the long range tail. What is 
important is that at the size of the barrier |a: — j/| ^ I <Si L the interaction ( |l8|) has the right form, V = ll\x — y\. 

In the HF method the many body wave function of the system (|l^) is represented in the form of Slater 
determinant of single particle orbitals ipi{x). Each orbital obeys the equation 
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h(pi = ei(pi, 



(21) 



where is the single particle energy and h is the HF Hamiltonian 



h(p{x) 



(22) 



Ueff^U{x)+J2 I \Vjiy)\^V{x,y)dy. 



Here the summation is performed over all filled orbitals. To avoid misunderstanding we have to comment on 
a terminological question. Usually only the case of zero vector potential, ^ = 0, is called the HF method. 
li A the method is called the Hartree-Fock method in the external field, or if A is time dependent it is 
called the time dependent Hartree-Fock methodE2l. The external field method is equivalent to the Random 
Phase Approximation (RPA). A crucial point is that the external field Hartree-Fock metiiod (=RPA) is gauge 
invariant, and therefore the electric current calculated within this method is conservedEj. 

For computations we use a ring of length L = 80 with a finite grid of 400 or 600 points. A naive finite 
grid version of the Hamiltonian (|2^) has a weak dependence on the gauge and hence the electric current is 
not exactly conserved. To avoid this minor trouble we use the lattice (=grid) modification of the Hamiltonian 
(H) replacing (p - A)'^ip by [2|<y9(n)p - ip*{n + l)e*"^''(p(n) - ip* {n)e~''^'^Lp{n + 1)] /2/i^. Here h is the step of 
the grid and ip{n) is the wave function on the site n of the grid. The electric current corresponding to the grid 
Hamiltonian is 



J = 



E 4 We^-'V.C" + 1) - ^*in + l)e-^V. W] • 



(23) 



This current is exactly conserved at a finite h, and at ft, it coincides with the standard electric current. 

First we consider the impurity pinned Wigner crystallizationEj. The electron density on the ring found as 
a result of the self consistent solution of the HF equations without any impurity, U{x) = 0, is plotted in Fig.l 



by long dashed lines. We consider the total number of electrons 
density is homogeneous. The density with the external potential 



U{x) = 







g(|x|-i/2)/d _^ I 

at C/q = 0.5, / — 8, and d—\ are shown by solid lines. 



N^+N^^ 78, 50, 38, 18, A^t = ^i- The 



(24) 
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FIG. 1. The Hartree-Fock electron density on the ring without (long dashed line) and with (solid line) small impurity. 
The dashed line corresponds to the Hartree approximation. The number of electrons on the ring is N — 78, 50, 38, 18. 
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It is clearly seen that the impurity pinned charge density wave is developed at the electron density n « 0.5. 
For comparison, we also show by dashed lines the results obtained in the Hartree approximation (no exchange 
interaction). The Hartree approximation also gives the charge density wave but underestimates the effect. 
In real Wigner crystal spins are localized at the sites of the lattice. This spin structure certainly can not be 
obtained in the HF approximation that enforces zero spin density at any point. However we are not going to 
study a real Wigner crystal at very low electron density. Anyway such a state does not conduct an electric 
current. We are interested in n > 0.5. The charge density wave at n ~ 0.5 is only a precursor for the Wigner 
crystallization and this precursor can be described in the HF approximation. 

In obtaining the results presented on Fig.l we have used the transverse cutoff at = 2, see eqs. (p^),(pO|). 
At the smaller value of at the charge density wave is developed at a higher electron density. This reflects the 
fact that the HF method with the ideal Coulomb interaction overestimates the tendency towards the charge 
density waves. Fottunately the value of at ~ 2 compensates this overestimation giving the correct value of the 
"critical" densityllj. This is why the application of the Hartee-Fock method is justified. 

We also performed HF calculations with total spin S ^ 0. Technically it is very simple, we just set ^ -/Vj. 
The macroscopic magnetization, i.e. S' ~ TV, does not arise. However, it is interesting to note that at the low 
electron density, n < 1, the ground state with = 1 or S" = 2 sometimes has slijglitly lower energy than the 
state with S ^ 0. This clearly contradicts to the rigorous Lieb-Mattis theorerrO and therefore it is just a 
byproduct of the HF approximation. We disregard this result and consider only the states with S = 0. 



IV. MANY-BODY CALCULATION OF THE BARRIER TRANSMISSION PROBABILITY 



For calculations in this section we take the number of electrons on the ring N ^ N]- + Ni = 158, = A^i- 
This corresponds to the average electron density < n >= 1.975 that is well above the charge density wave 
threshold. The extfijrnal potential barrier is taken in the form ( B4|). The parameter Uq models the gate potential 
in the experimentaJ^lZl. Selfconsistent solution of the HF eqs. (^1|) is performed at J7o = and at f/o ^ and 
then the transmission probability is found using eq. (^6|) . The transmission probability as a function of Uq at 
d — 0.5 and the values of the barrier length i = 4, 6, 8, 10, 12 is plotted in Fig. 2. 




FIG. 2. The transmission probability T versus the potential Uo- The transmission probability has been calculated in 
the Hartree-Fock approximation at d = 0.5 for the values of the barrier length 1 — 4,6, 8, 10, 12. 

The plot of T versus Uq is representative if L ^ /. However in the computations L = 80 and therefore the 
inequality is not so strong. As a result the Fermi energy attains a weak Uq dependence. To compensate this 
effect, in Fig. 3 we plot the transmission probability T versus the value ofVb — Uq — ep{Uo) + eF{Ua ~ 0) that 
models the gate potential more accurately. 
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FIG. 3. The transmission probability T versus the gate potential Vq. The transmission probability has been calculated 
in the Hartree-Fock approximation at d — 0.5 for the values of the barrier length 1 — 4,6, 8, 10, 12. 



There is no a qualitative difference between Fig. 2 and Fig. 3, but Fig. 3 is more correct quantitatively. There 
are some very sharp structures in Fig. 2 and Fig. 3, but they arise only because of the relatively large step in 
Uq. The computations are rather time consuming and this limits the number of points. For comparison we 
present in Fig. 4 the same plots as in Fig. 3, but the curves have been obtained in the Hartree approximation. 
So the exchange term in eq. (^) has been dropped out. 
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FIG. 4. The transmission probability T versus the gate potential Vo- The transmission probability has been calculated 
in the Hartree approximation at d — 0.5 for the values of the barrier length 1 — 4,6, 8, 10, 12. 

The results of Hartree approximation are qualitatively similar to that of the HF approximation, Fig. 3. However 
all the structures are less pronounced. 

The previous figures present the barrier transmission probability calculated at the parameter d = 0.5, see eq. 
(p^). In Fig. 5 and Fig. 6 we plot the similar results at d = 1. The Fig. 5 corresponds to the HF approximation 
and Fig. 6 corresponds to the Hartree approximation. 
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FIG. 5. The transmission probability T versus the gate potential Vq. The transmission probability has been calculated 
in the Hartree-Fock approximation at d — 1 for the values of the barrier length I — 8, 10, 12, 14, 16. 
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FIG. 6. The transmission probability T versus the gate potential Vo- The transmission probability has been calculated 
in the Hartree approximation at d — 1 for the values of the barrier length I = 8, 10, 12, 14, 16. 



The Hartree approximation gives results qualitatively similar to that of the HF approximation, but the Hartree 
approximation clearly underestimates the structure. 

Fig. 3 and Fig. 5 clearly demonstrate a shoulder at the value of the transmission probability T ^ 0.7. For 
longer Jparriers an additional structure appears at smaller values of T. This is in agreement with experimental 
datcoill. We remind that we use the unit of length as ~ 10^^/im, so the typical length of barriers presented 
in Fig. 3 and Fig. 5 is between O.l/im and 0.2/im. 

To understand the reason for the structures in the transmission probability we plot the effective selfconsistent 
potential C4// (see eq. (22)) and the electron density above the barrier at parameters corresponding to these 
structures. The point A in Fig. 5 indicates the shoulder at T « 0.65 for the barrier of the length I = 10, the 
parameters are: Z = 10, d = 1, Uq — 9.8. The corresponding selfconsistent potential is shown in Fig. 7 by the 
solid line, the external potential (|2j) is shown by the dashed line. 
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FIG. 7. The selfconsistent HF potential \23j is shown by the solid line. The external potential 
dashed line. The parameters corresponds to the point A in Fig. 5: I = 10, d = 1, Uo = 9.8. 



is shown by the 



For an estimate the external potential in Fig. 7 can be approximated near the top by a parabola, U ■ 
0.5. 



with uj 

a\\/at ^ ^Ju}±/uj\\ 



This gives the ratio of the typical longitudinal size to the typical transverse size of the contact, 
1, see also eq. (|20|). However the selfconsistent potential. Fig. 7, is very much different 
from the simple parabola because of redistribution of electrons along the contact. The selfconsistent potential 
calculated in the Hartree approximation is practically the same as that calculated in the HF approximation. 
Nevertheless the HF transmission probability has a shoulder at the point A, see Fig. 5, and the Hartree 
transmission probability has no a shoulder, see Fig. 6. The difference is due to the exchange interaction. 
To demonstrate this we plot in Fig. 8 the electron density above the barrier. The HF electron density is shown 
by the solid line, and the Hartree electron density is shown by the dashed line. We see that the shoulder in the 
transmission probability is related to the charge density wave above the barrier. There is no a charge density 
wave in the Hartree approximation, and there is no a shoulder in the corresponding transmission probability 
shown in Fig. 6. 
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FIG. 8. The electron density on the barrier at I = 10, d = 1, Uo = 9.8. This corresponds to the shoulder A in 
Fig. 5. The solid line is obtained in the Hartree-Fock approximation and the dashed line is obtained in the Hartree 
approximation. 
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The point C in Fig. 5 indicates the local minimum at T « 0.5 for the barrier of the length I = 16, the 
parameters are: I = 16, d = 1, Uq = 11.4. The corresponding selfconsistent potential is shown in Fig. 9 by the 
solid line, the external potential pi) is shown by the dashed line. 
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FIG. 9. The selfconsistent HF potential i\23[ ) is shown by the solid line. The external potential 
dashed line. The parameters corresponds to the point C in Fig. 5: I — 16, d = 1, i7o = 11.4. 



is shown by the 



Similar to the previous case the HF selfconsistent potential is practically indistinguishable from the Hartree 
one. Nevertheless the HF transmission probability has a well pronounced minimum, see point C in Fig. 5, and 
the Hartree transmission probability has just a shoulder, see Fig. 6. The difference is due to the exchange 
interaction. The HF electron density is shown in Fig. 10 by the solid line, and the Hartree electron density is 
shown by the dashed line. There is an additional modulation of the HF density due to the exchange interaction. 
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FIG. 10. The electron density on the barrier at I = 16, d = 1, Uo = 11.4. These parameters correspond to the local 
minimum C in Fig. 5. The solid line corresponds to the Hartree-Fock approximation and the dashed line corresponds to 
the Hartree approximation. 



We see that the structures in the transmission probability are related to the resonance structures above the 
potential barrier. The resonance structures are related to the charge density waves developing on the barrier 
when the average electron density above the barrier is of the order of O.S/as. 
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V. CONCLUSIONS 



The fictitious gauge field metliod lias been developed to calculate a potential barrier transmission probability 
T with account of many-body effects. The transmission probability is directly related to the conductance, 
G = The results of the Hartree-Fock calculations of the transmission probability are shown in Fig. 3 and 

Fig. 5. The results demonstrate a plateau or a shoulder at the value of the transmission probability T w 0.7. 
For longer barriers this structure is getting weak, but additional structures appear at lower values of T . All 
the structures are related to the development of the charge density waves on the barrier. This is a precursor 
for Wigner crystallizatie]*. We believe that this explains the conductance anomalies observed in experiments 
with quantum contactsai3. 

We would like to note that there are two further questions that can be studied within the developed fictitious 
gauge field approach. The first question is how the structures in the conductance evolve in the external magnetic 
field. The experimental dataa indicates that the conducatnce anomaly is somehow related to the spin. In the 
picture considered in the present work the effect is certainly spin dependent because it is due to the exchange 
interaction. To analyze this problem in details one should perform calculations at 7^ N^, so the problem is 
more technically involved compared to the considered case. Another question is the temperature dependence 
of the anomalies. It is known from experiment that the effect is increasing with temperature up to ~ IK. 
Unfortunately there is no a direct way to consider the finite temperature case within the developed formalism. 
However it is clear that any temperature dependence is related to the excitations above the barrier. Moreover 
there is a direct spectroscopic observation of such excitationQ. The problem of the collective excitation can 
be addresses within the developed formalism. To approach this problem one should apply the time dependent 
Hartree-Fock method in the external gauge field. This calculation is substantially more technically involved 
compared to the considered case. 

I am grateful to P. E. Lindelof who has attracted my attention to this problem. I am also grateful to V. 
V. Flambaum, C. J. Hamer, M. I. Katsnelson, M. Yu. Kuchiev, A. Nersesyan, D. J. Reilly, and J. Thakur for 
helpful conversations. 
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